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Abstract — A single source network is said to be memory-free if 
all of the internal nodes (those except the source and the sinks) 
do not employ memory but merely send linear combinations of 
the symbols received at their incoming edges on their outgoing 
edges. In this work, we introduce network-error correction for 
single source, acyclic, unit-delay, memory-free networks with 
coherent network coding for multicast. A convolutional code is 
designed at the source based on the network code in order to 
correct network-errors that correspond to any of a given set 
of error patterns, as long as consecutive errors are separated 
by a certain interval which depends on the convolutional code 
selected. Bounds on this interval and the held size required 
for constructing the convolutional code with the required free 
distance are also obtained. We illustrate the performance of 
convolutional network error correcting codes (CNECCs) designed 
for the unit-delay networks using simulations of CNECCs on an 
example network under a probabilistic error model. 

I. Introduction 

Network coding was introduced in [1] as a means to 
improve the rate of transmission in networks, and often achieve 
capacity in the case of single source networks. Linear network 
coding was introduced in [2]. An algebraic formulation of 
network coding was discussed in [3] for both instantaneous 
networks and networks with delays. 

Network error correction, which involved a trade-off be- 
tween the rate of transmission and the number of correctable 
network-edge errors, was introduced in [4] as an extension 
of classical error correction to a general network setting. 
Along with subsequent works [5] and [6], this generalized the 
classical notions of the Hamming weight, Hamming distance, 
minimum distance and various classical error control coding 
bounds to their network counterparts. In all of these works, it 
is assumed that the sinks and the source know the network 
topology and the network code, which is referred to as 
coherent network coding. Network error correcting codes were 
also developed for non-coherent (channel oblivious) network 
coding in [7], [8] and [9]. Network error correction under 
probabilistic error settings has been studied in [10]. Most 
recently, multishot subspace codes were introduced in [11] for 
the subspace channel [7] based on block-coded modulation. 

A set of code symbols generated at the source at any 
particular time instant is called a generation of code symbols. 
So far, network error correcting schemes have been studied 
only for acyclic instantaneous (delay-free) networks in which 
each node could take a linear combination of symbols of only 
the same generation. 

Convolutional network codes were discussed in [12]— [14] 
and a connection between network coding and convolutional 
coding was analyzed in [15]. Convolutional network error 



correcting codes (which we shall henceforth refer to as 
CNECCs) have been employed for network error correction 
in instantaneous networks in [16]. 

A network use [16] is a single usage of all the edges of 
the network to multicast utmost min-cut number of symbols 
to each of the sinks. An error pattern is a subset of the set of 
edges of the network which are in error. It was shown in [16] 
that any network error which has its error pattern amongst a 
given set of error patterns can be corrected by a proper choice 
of a convolutional code at the source, as long as consecutive 
network errors are separated by a certain number of network 
uses. Bounds were derived on the field size for the construction 
of such CNECCs, and on the minimum separation in network 
uses required between any two network errors for them to be 
correctable. 

Unit-delay networks [13] are those in which every link 
between two nodes has a single unit of delay associated with 
it. In this work, we generalize the approach of [16] to the case 
of network error correction for acyclic, unit-delay, memory- 
free networks. We consider single source acyclic, unit-delay, 
memory-free networks where coherent network coding (for the 
purpose of multicasting information to a set of sinks) has been 
implemented and thereby address the following problem. 

Given an acyclic, unit-delay, single source, memory-free 
network with a linear multicast network code, and a set of 
error patterns $, how to design a convolutional code at the 
source which will correct network errors corresponding to 
the error patterns in $, as long as consecutive errors are 
separated by a certain number of network uses? 

The main contributions of this paper are as follows. 

• Network error correcting codes for unit-delay, memory- 
free networks are discussed for the first time. 

• A convolutional code construction for the given acyclic, 
unit-delay, memory-free network that corrects a given 
pattern of network errors (provided that the occurrence 
of consecutive errors is separated by certain number of 
network uses) is given. For the same network, if the 
network code is changed, then the convolutional code 
obtained through our construction algorithm may also 
change. Several results of this paper can be treated as 
a generalization of those in [16]. 

• We derive a bound on the minimum field size required 
for the construction of CNECCs for unit-delay networks 
with the required minimum distance, following a similar 
approach as in [16]. 

• We also derive a bound on the minimum number of 
network uses that two error events must be separated by 
in order that they get corrected. 
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> We also introduce processing functions at the sinks in 
order to address the realizability issues that arise in the 
decoding of CNECCs for unit-delay networks. 

• We show that the unit-delay network demands a CNECC 
whose free distance should be at least as much as that of 
any CNECC for the corresponding instantaneous network 
to correct the same number of network errors. 

• Using a probabilistic error model on a modified butterfly 
unit-delay memory-free network, we use simulations to 
study the performance of different CNECCs. 

• Towards achieving convolutional network error correc- 
tion, we address the issue of network coding for an 
acyclic, unit-delay, memory-free network. As a by- 
product, we prove that an n-dimensional linear network 
code (a set of local kernels at the nodes) for an acyclic, 
instantaneous network continues to be an n-dimensional 
linear network code (i.e the dimension does not reduce) 
for the same acyclic network, however being of unit-delay 
and memory-free nature. 

The rest of the paper is organized as follows. In Section HH 
we discuss the general network coding set-up and network 
errors. In Section [TTTJ we give a construction for an input 
convolutional code for the given acyclic, unit-delay, memory- 
free network which shall correct errors corresponding to a 
given set of error patterns and also derive some bounds on the 
field size and minimum separation in network uses between 
two correctable network errors. In Section IIV1 we give some 
examples for this construction. In Section [V] we provide 
a comparison between CNECCs for instantaneous networks 
[16] and those for unit-delay, memory-free networks of this 
paper. In Section IVI1 we discuss the results of simulations 
of different CNECCs run on a modified butterfly network 
assuming a probabilistic model on edge errors in the network. 
We conclude this paper in Section fVTll with some remarks and 
some directions for further research. 

II. Problem Formulation - CNECCs for unit-delay, 

MEMORY-FREE NETWORKS 

A. Network model 

We consider acyclic networks with delays in this paper, the 
model for which is as in [3], [13]. An acyclic network can be 
represented as an acyclic directed multi-graph (a graph that 
can have parallel edges between nodes) Q = (V, £) where V 
is the set of all vertices and £ is the set of all edges in the 
network. 

We assume that every edge in the directed multi-graph 
representing the network has unit capacity (can carry utmost 
one symbol from ¥ q ). Network links with capacities greater 
than unit are modeled as parallel edges. The network has 
delays, i.e, every edge in the directed graph representing the 
input has a unit delay associated with it, represented by the 
parameter z. Such networks are known as unit-delay networks. 
Those network links with delays greater than unit are modeled 
as serially concatenated edges in the directed multi-graph. The 
nodes of the network may receive information of different 
generations on their incoming edges at every time instant. We 
assume that the internal nodes are memory-free and merely 



transmit a linear combination of the incoming symbols on their 
outgoing edges. 

Let s € V be the source node and T be the set of all 
receivers. Let n T be the unicast capacity for a sink node TeT 
i.e the maximum number of edge-disjoint paths from s to T. 
Then 

n = min n„ 
TeT 

is the max-flow min-cut capacity of the multicast connection. 

B. Network code 

We follow [3] in describing the network code. For each node 
v E V, let the set of all incoming edges be denoted by Ti(v). 
Then |Tj (v) | = 6t(v) is the in-degree of v. Similarly the set 
of all outgoing edges is defined by To(v), and the out-degree 
of the node v is given by |To(v)| = 5o{v). For any e G £ and 
v G V, let head(e) = v, if v is such that e G Similarly, 
let tail(e) = v, if v is such that e G To(v). We will assume 
an ancestral ordering on £ of the acyclic graph Q. 

The network code can be defined by the local kernel 
matrices of size 8i(v) x So(v) for each node v G V with 
entries from ¥ q . The global encoding kernels for each edge 
can be recursively calculated from these local kernels. 

The network transfer matrix, which governs the input-output 
relationship in the network, is defined as given in [3] for an 
n dimensional network code. Towards this end, the matrices 
A,K,and _B T (for every sink T G T are defined as follows. 

The entries of the n x \£ | matrix A are defined as 

A = f oti jej if e 3 G r (s) 
1,3 1 otherwise 

where a,- )e , G ¥ q is the local encoding kernel coefficient at 
the source coupling input i with edge ej G To(s). 

The entries of the \£\ x \£\ matrix K are defined as 

_ ( (3i,j if head(ei) = tail(ej) 
l ' 3 1 otherwise 

where the set of fli.j G ¥ q is the local encoding 
kernel coefficient between e; and ej, at the node 

v = head(ei) — tail(ej). 

For every sink T G T, the entries of the \£\ x n matrix B T 
are defined as 

/,,/ f >■ ,i if ejtTifT) 
l >i \ otherwise 

where all e e ,i G ¥ q . 

For unit-delay, memory-free networks, we have 

F(z) := (I -zKy 1 

where / is the \£\ x \£\ identity matrix. Now we have the 
following definition. 

Definition 1 ( [3]): The network transfer matrix, Mt{z), 
corresponding to a sink node T G T is a full rank (over 
F 9 (z)) n x n matrix defined as 

M T (z) := AF(z)B T = AF T (z). 
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With an n-dimensional network code, the input and the 
output of the network are n-tuples of elements from F 9 [[z]]. 
Definition [T] implies that if x(z) G [[2;]] is the input to 
the unit-delay, memory-free network, then at any particular 
sink T G T, we have the output, y(z) G F"[[z]], to be 
y{z)=x{z)M T {z). 

C. CNECCs for single source, unit-delay, memory-free net- 
works 

A primer on the basics of convolutional codes can be 
found in Appendix [A] Assuming that an n-dimensional linear 
network code multicast has been implemented in the given 
single source unit-delay, memory-free network, we extend the 
definitions of the input and output convolutional codes of 
CNECCs for instantaneous networks from [16] to the unit- 
delay, memory-free case. 

Definition 2: An input convolutional code, C s , correspond- 
ing to an acyclic, unit-delay, memory-free network is a con- 
volutional code of rate k/n(k < n) with a input generator 
matrix Gi(z) implemented at the source of the network. 

Definition 3: The output convolutional code Ct, corre- 
sponding to a sink node T G T in the acyclic, unit-delay, 
memory-free network is the k/n(k < n) convolutional code 
generated by the output generator matrix Go.t(z) which is 
given by 

Go,t(z) = G i (z)Mt(z) 

with Mt(z) being the full rank network transfer matrix 
corresponding to an n-dimensional network code. 




Fig. 1. A simple min-cut 2 network with one source and one sink 

Example 1: Consider the single source, single sink network 
as shown in FigQ] Let the field under consideration be Fa. The 
local kernels at the intermediate node are unity. Therefore the 
network transfer matrix at the sink is (assuming the given 
ancestral ordering) 



M T (z) = 



z 
z 2 



Suppose we choose the input convolutional code C s to be 
generated by the matrix 



Gt(z) 



\l + z 2 1 + z + z 2 ] 



Then the output convolutional code Ct is generated by 

Go,t(z) = [z + z 3 z 2 + z 3 + z 4 ] . 

D. Network errors 

Observing a 'snap-shot' of the network at any particular 
time instant, we define the following terms. An error pattern 
p, as stated previously, is a subset of £ which indicates the 
edges of the network in error. An error vector w is a 1 X \£ \ 
vector which indicates the error occurred at each edge. An 
error vector is said to match an error pattern (i.e w G p) if 
all non-zero components of w occur only on the edges in p. 
An error pattern set $ is a collection of subsets of £ , each of 
which is an error pattern. 

Let x{z) G F™[[z]] be the input to the network , and w G 

\£\ 

Fq 1 be the error vector corresponding to the network errors 
that occurred at any time instant i (i G Zjj\ referenced from 
the first input time instant). Then, the output, y(z) G F"[[z]] 
at any particular sink T G T can be expressed as 

y(z) = x(z)M T (z) + wz i F T (z). 

In case there are a number of errors at a number of time 
instants, we have the formulation as 

y(z) = x(z)M T (z) + w{z)F T (z) 

I £ 

wherein every monomial of w(z) G FL 1 [[z]] of the form WiZ 1 
incorporates the error vector Wi G F q occurring at the time 
instant i. 

iii. cneccs for unit-delay, memory-free 
networks - Code Construction and Capability 

A. Network code for acyclic unit-delay memory-free networks 

In Section IIII-BI we give a construction of a CNECC for a 
given acyclic, unit-delay, memory-free network. Towards that 
end, we first address the problem of constructing network 
codes for acyclic, unit-delay, memory-free networks. Although 
network code constructions have been given for acyclic instan- 
taneous networks [17], the problem of constructing network 
codes for acyclic, unit-delay, memory-free networks is not 
directly addressed. The following lemma shows that solving 
an n-dimensional network code design problem for an acyclic, 
unit-delay, memory-free network is equivalent to solving that 
of the corresponding acyclic instantaneous network with the 
same number of dimensions. 

Lemma 1: Let Q(V,£) be a single source acyclic, unit- 
delay, memory-free network, and Ginst be the corresponding 
instantaneous network (i.e with the same graph as that of Q, 
but no delay associated with the edges). Let Af be the set of all 
5 1 (v) x 5o(v) matrices V v G V, i.e, the set of local encoding 
kernel matrices at each node, describing an n'-dimensional 
network code (over F 9 ) for Qi ns t in 1 < min-cut of the source- 
sink connections in Ginst)- Then the network code described 
by Af continues to be an n'-dimensional network code (over 
F g (z)) for the unit-delay, memory-free network Q. 

Proof: Let Mt be the n' x n' network transfer matrix of 
any particular sink node T G T in Qi ns t, and Mt{z) be the 
n! x n' network transfer matrix of the same sink T in Q. We 
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first note that the matrix Mt can be obtained from Mt(z) by 
substituting z = z° = 1, i.e, 

M T = M T {z)\ z=1 . 

Given that Mt is full rank over ¥ q , we will prove that Mt(z) 
is full rank over ¥ q (z) by contradiction. 

Suppose that Mt(z) was not full rank over ¥ q (z), then we 
will have 

i—n —1 /\ 

E -rTT m i( z ) = m n>(z) (1) 
t=l l ^ z ' 

where mi(z) is the i th row of Mt{z) and a.i(z), bi(z) G F g [z] 

V i = 1, 2, .., n' are such that bi(z) ^ 0, a.j(z) 7^ for at least 
one i, and gcd(a,i(z) , bi(z)) = 1, V i. 

We have the following two cases 
Case 1: bi{z)\ z —\ ^ Vi. 
Substituting z — 1 in ([T), we have 

z— n ; — 1 

E y^rnj = m„/ (2) 

where a* = a t (z)\ z=1 , b { = bi(z)\ z=1 and rrii = mi(z)\ z=1 
is the i th row of Mr. 

Clearly m n i ^ since My is full rank, and hence the left 
hand side of (f2]i can't be zero. Therefore some non-zero linear 

th 

combination of the first n' — 1 rows of Mt is equal to its n' 
row, which contradicts the given statement that Mt is full 
rank over ¥ q . Therefore Mr (2) must be full rank over ¥ q (z). 

Case 2: bi(z)\ z =i = for at least one i. 

Let T C {1,2, ...,n'} such that (z - l) p '\bi(z) for some 
positive integer p' . Let p be an integer such that 

p = maxp' 
iei' 

Now, from (Q~|i we haven 

E (z-l) p ^rlm l (z) = (z-irm n/ (z) (3) 

Let I C {1,2, ..,71'} such that (z - l) p |foi(z) Vi E 1 
Then we must have that (2 — 1) { aj(z) V i € X, since 
gcd( ai {z), h(z)) = 1. Also, let t;(z)^,(z)/( z -lfeF,[z] 

V i £ I. Hence we have 

(<-<g)^-(§$)l~-^'.\<°>- 

where b\ = b[{z)\ z=l e ¥ q \{0}, since (z - 1) f &J(jb). 
Substituting z = 1 in (|3), we have 

E = 

i.e, a non-zero linear combination of the rows of My is 
equal to zero, which contradicts the full-rankness of Mt, thus 
proving that Mt(z) has to be full rank over ¥ q (z). ■ 



B. Construction 

This subsection presents the main contribution of this work. 
We assume an n dimensional network code (n being the min- 
cut) on this network has implemented on the given network 
which is used to multicast information to a set of sinks. We 
describe a construction of an input convolutional code for 
the given acyclic, unit-delay, memory-free network which can 
correct network errors with patterns in a given error pattern 
set, as long as they are separated by certain number of network 
uses. 

Let Mt(z) = AFt(z) be the nxn network transfer matrix 
from the source to any particular sink T G T. Let $ be the 
error pattern set given. We then define the processing matrix 
at sink T, Pt(z), to be a polynomial matrix as 

Pt{z)^p t {z)M^{z) 

where p T (z) G ¥ q [z] is some processing function chosen 
such that Pt{z) is a polynomial matrix. Now, we have the 
construction of a CNECC for the given network as follows. 

1 ) We first compute the set of all error vectors having their 
error pattern in $ that is defined as follows 

W * = U { w = (wi,W2,-,w\ £ \) e Fjf \wep\. 

2) Let 

W T ■= {wF T (z) I w G W.J,} (4) 

be computed for each sink T. This is the set of rt-tuples 
(with elements from F g [z]) at the sink T due to errors 
in the given error patterns $. 

3) Let the set W s C F"[z] 

W s := |J {w T (z)P T (z) I w T (z) G W T } (5) 
TeT 

be computed. 

4) Let 

t s = max wh (wJz)) . 

where wh indicates the Hamming weight over ¥ q . 

5) Choose an input convolutional code C s with free distance 
at least 2t s + 1 as the CNECC for the given network. 

C. Decoding 

Before we discuss the decoding of CNECCs designed ac- 
cording to Subsection IIII-BI we state some of the results from 
[16] related to the bounded distance decoding of convolutional 
codes in this section. 

Let C be a rate b/c convolutional code with a gener- 
ator matrix G(z). Then, corresponding to the information 
sequence tto, u\, ..(ttj G F*) and the codeword sequence 
Vo,Vi, ...(v.i G Fg), we can associate an encoder state se- 
quence cr , <xi, .., where <r t indicates the content of the delay 
elements in the encoder at a time t. We define the set of j 
output symbols as 

V[0,j) •= [V0,V1, ...,Vj-i] 
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The parameter Td 



(C) [16] is defined as follows. 

i + i 



Td free (C) := max 

where Sd frcc [16] is defined as the set of all possible truncated 
codeword sequences f[o,j) of weight less than df ree (C) that 
start in the zero state is defined as follows 

Sd free ■= {v[oj) | w H (v[ j)) < df ree (C),(To = 0,V j > 0} 

where iwjj indicates the Hamming weight over ¥ q . 

Then, we have the following proposition. 

Proposition 1 ( [16] ): The minimum Hamming weight 
trellis decoding algorithm can correct all error sequences 
which have the property that the Hamming weight of the 
error sequence in any consecutive Td free (C) segments (a 
segment being the set of c code symbols generated for every 
c information symbols) is utmost d f r 'e(fi) — ! 

Now, we discuss the decoding of CNECCs for unit-delay 
memory-free networks. Let Gj(z) be the k X n generator 
matrix of the input convolutional code, C s , obtained from 
the given construction. Let Go,t(z) = G i(z)Mt(z) be the 
generator matrix of the output convolutional code, Ct, at sink 
T G T, with Mt(z) being its network transfer matrix. 

For each sink TeT, let 

tx = max wh{w t {z)). 

w t (z)£Wt 

Let rriT be the largest integer such that 

dfree(C T ) > 2m T t T + 1. (6) 

Clearly, rax > 0. Each sink can choose decoding on the trellis 
of the input or its output convolutional code based on the 
characteristics of the output convolutional code as follows 

Case-A: This is applicable in the event of all of the 
following conditions being satisfied. 

i.) 

m T > 1 (7) 



ii.) 



7d /ree (Cr) < rriTTd fre<s (C s ). 



(8) 



iii.) The output convolutional code generator matrix 
Go,t(z) is non-catastrophic. /m 

In this case, the sink T performs minimum distance decoding 
directly on the trellis of the output convolutional code, Ct- 

Case-B: This is applicable if at least one of the 3 conditions 
of Case-A is not satisfied, i.e, if either of the following 
conditions hold 

i. ) m,T = 

ii. ) m T > 1 and T dfree (C T ) > mrTd freB (C s ). 

iii. ) The output convolutional code generator matrix 

Go,t{z) is catastrophic. 
This method involves processing (matrix multiplication using 
Pt{z)) at the sink T. We have the following formulation at 
the sink T. Let 



K( z ) v' 2 (z) ... v' n (z)] 



[vi(z) v 2 (z) ... v n (z)} 
+ K(z) w 2 (z) ... w n (z)} 



represent the output sequences at sink T, where 

[vi(z) v 2 (z) ... v n (z)] 

= u(z)G ,t(z) = u(z)G I (z)M T (z) 

u{z) being the k length vector of input sequences, and 

[wi(z) w 2 (z) ... w n (z)} 

represent the corresponding error sequences. Now, the output 
sequences are multiplied with the processing of the network 
transfer matrix Pt(z), so that decoding can be done on the 
trellis of the input convolutional code. Hence, we have 

[<(*) v' 2 \z) ... <(*)] 

= Wii z ) v' 2 {z) ... v' n {z)]P T {z) 

= u(z)p T (z)G I {z) + [w 1 (z) w 2 (z) ... w n (z)]P T (z) 

= u(z)p T (z)G I (z) + [w[(z) w' 2 (z) ... w' n (z)} 

where w'(z) = [w'^z) w' 2 (z) ... w' n (z)] now indicate the 
set of modified error sequences that are to be corrected. Now 
the sink T decodes to the minimum distance path on the 
trellis of the code generated by p T (z)Gi(z), which is the input 
convolutional code as Gi(z) and p T (z)Gi(z) are equivalent 
generator matrices. 

Remark 1: In [16], the approach to the construction of a 
CNECC for an instantaneous network was the same as in here. 
However, the set W s was defined in [16] as 



W s := |J {w T M~ 1 | w T e W T } 



= [j {wFtM- 1 I wep} (11) 

where the network transfer matrix Mt and Ft correspond to 
a sink T in the instantaneous network. 

In this paper, the definition for W s is as in (0 and involves 
the processing matrix Pr(z) instead of the inverse of the net- 
work transfer matrix. The processing function p T (z) for a sink 
T is introduced because of the fact that the matrix M^ 1 (z) 
might not be realizable and also for easily obtaining the 
Hamming weight of the error vector reflections (w s (z) G W s ) 
by removing rational functions in M^ 1 (z). 

The degree of the processing function p T (z) directly in- 
fluences the memory requirements at the sinks and therefore 
should be kept as minimal as possible. Therefore, with 



M T \z) 



M(z) 



Det (M T (z)) 



where the n x n matrix M(z) is the adjoint of Mt(z), ideally 
we may choose p T (z) as follows. 



p T {z) = 



Det {M T (z)) 
9{z) 



(12) 



where g(z) = gcd(rriij(z), V 1 < j,j < n), rriij(z) being 
the (i,j) th element of M(z). 
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D. Error correcting capability 

In this subsection we prove a main result of the paper 
given by Theorem Q] which characterizes the error correcting 
capability of the code obtained via the construction of Subsec- 
tion IIII-BI We recall the following observation that in every 
network use, n encoded symbols which is equal to the number 
of symbols corresponding to one segment of the trellis, are to 
be multicast to the sinks. 

Theorem 1: The code C s resulting from the construction of 
Subsection IIII-BI can correct all network errors that have their 
pattern as some p £ $ as long as any two consecutive network 
errors are separated by T^, ree (C s ) network uses. 

Proof: We first prove the theorem in the event of Case- 
A of the decoding. Suppose the network errors are such 
that consecutive network errors are separated by Td /rec (C s ) 
network uses. Then the vector of error sequences at sink T, 
w T (z), is such that in every T^, ree (C s ) segments, the error 
sequence has utmost tr Hamming weight (over F g ). Therefore 
in rriTTd fm (C s ) segments, the Hamming weight of the error 
sequence would be utmost m^T- 

Then the given condition ((HJ would imply that in every 
Td fm (Ct) segments of the output trellis, the error sequences 
have Hamming weight utmost tot^t- Condition (O together 
with © and Proposition Q] implies that these error sequences 
are correctable. This proves the given claim that errors with 
their error pattern in $ will be corrected as long as no two 
consecutive error events occur within Td }rcc (C s ) network uses. 

In fact, condition © and © implies that network errors 
with pattern in $ will be corrected at sink T, as long as 
consecutive error events are separated by T^. e (Cr)- 

Now we consider Case B of the decoding. Suppose that the 
set of error sequences in the formulation given, w'(z), is due 
to network errors that have their pattern as some such 
that any two consecutive such network errors are separated by 
at least 7d /rcc (C s ) network uses. 

Therefore, along with step 4 of the construction, we have 
that the maximum Hamming weight of the error sequence 
w'{z) in any consecutive Td frcc (C s ) segments (network uses) 
would be utmost t s . Because of the free distance of the code 
chosen and along with Proposition Q] we have that such errors 
will get corrected when decoding on the trellis of the input 
convolutional code. ■ 

E. Bounds on the field size and Id /ree (C s ) 

1 ) Bound on field size: Towards obtaining a bound on the 
sufficient field size for the construction of a CNECC meeting 
our free distance requirement, we first prove the following 
lemmas. 

Lemma 2: Given an acyclic, unit-delay, memory-free net- 
work Q(V, £) with a given error pattern set let Tdeiay — 1 be 
the maximum degree of any polynomial in the F(z) matrix. 
Let w h indicate the Hamming weight over ¥ q . If r is the 
maximum number of non-zero coefficients of the polynomials 
p T (z) corresponding to all sinks in T, i.e 

r = maxuifl- (p T (z)) , 



then 

max w H (w s (z)) < rn [{n + 1) {Tdeiay - 1) + 1] . 
w s (z)ew s 

where W s is as in © in Subsection IIII-BI 

Proof: Any element w s (z) £ ~W S indicates the n length 
sequences that would result in an output vector w T (z) at some 
sink T as a result of an error vector w in the network at time 
0, i.e 

w s (z) = wFt(z)p t (z)M^ 1 (z) = w T (z)p T (z)M^ 1 (z) 

Because of the fact that any polynomial in F(z) has degree 
utmost Tdeiay — 1> any error vector w at time can result in 
non-zero symbols (over F") in w T (z) at any sink T from the 
th time instant only upto utmost Tdeiay ~ 1 time instants. 

/T dclav -1 \ 

w T ( z ) = ( w t,i z J ■ 

where w T i £ F^ . 

The numerator polynomial of any element a(z) £ ¥ q (z) 
of the matrix M^ 1 {z) has degree utmost n {Tdeiay — 1). 
Therefore, considering the polynomial processing matrix 
Pt(z) = p T (z)M^ (z), we note that any element from 
Pt(z) has utmost r [n (Tdeiay — 1) + 1] non-zero components 
(over F g ), the worst case being r non-overlapping 'blocks' of 
n (Tdeiay — 1) + 1 non-zero components each. 

Therefore the first non-zero symbol of w T (z) (over F") at 
some time instant can result in utmost r [n (Tdeiay ~ 1) + 1] 
non-zero symbols in w s (z) (over F"). Henceforth, every 
consecutive non-zero symbol (over F") of w T (z) will result 
in utmost additional r F™ symbols in w s (z). Therefore any 
w 8 (z) £ W s is of the form 

/r[(n+l)(T dc!aB -l) + l] \ 

w s (z) = ^2 w Sti z % 

V i=o / 

where w s> i £ F™. Therefore the Hamming weight (over ¥ q ) 
of any w s (z) £ W s is utmost rn [(n + 1) (Tdeiay — 1) + 1], 
thus proving the lemma. ■ 
Our bound on the field size requirement of CNECCs for 
unit-delay networks is based on the bound on field size 
for the construction of Maximum Distance Separable (MDS) 
convolutional codes [18], a primer on which can be found in 
Appendix [B] 

Lemma 3: A (n, k) MDS convolutional code C (over some 
field Fg) with degree S — \(2t — l)k/n] can correct any 
error sequence which has the property that the Hamming 
weight(over ¥ q ) of the error sequence in any consecutive 
Td free (C) segments is utmost t. 

Proof: Because the generalized Singleton bound is satis- 
fied with equality by the MDS convolutional code, we have 

d fr ee(C) = (n - k) ( [5/k\ + 1 ) + 5 + 1 . 
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Substituting \(2t — 1) fc/n] for S, we have 

dfree (C) 

= (»-*) ( r(2 *~ fc 1)fc/n1 + 1) + r(« - 1) fc/»i + 1 



dfree{C) > (fl - k) ( 



(2t-l) 



1 



(2i- l)fc 



1 



d /ree (C) >2f + 1. 



Thus the free distance of the code C is at least It + 1, and 
therefore by Proposition [Tj such a code can correct all error 
sequences which have the property that in any consecutive 
Tdf rB e(C) segments, the Hamming weight (over ¥ q ) of the 
error sequence is utmost t. ■ 
For an MDS convolutional code being chosen as the input 
convolutional code (CNECC), we therefore have the following 
corollary 

Corollary 1: Let Q(V,£) be an acyclic, unit-delay, 
memory-free network with a network code over a suffi- 
ciently large field ¥ q and $ be an error pattern set, the 
errors corresponding to which are to be corrected. An (n, fc) 
input MDS convolutional code C s over ¥ q with degree 
5 = 2rk [(n + 1) (Tdeiay — 1) + 1] can be used to correct 
all network-errors with their error pattern in $ provided that 
consecutive network-errors are separated by at least Td free (C s ) 
network uses, where r and Tdeiay are as in Lemma [2] 

Proof: From Lemma [2] we have that in the construction 
of Subsection IIII-BI the maximum Hamming weight t s of any 
element in the set W s is utmost rn \(n + 1) (Tdeiay — 1) + 1] • 
For an input MDS convolutional code C s to be capable 
of correcting such errors with Hamming weight utmost 
rn [(n + 1) (Tdeiay — 1) + 1], according to Lemma|3] a degree 
5 = 2rk [(n + 1) (J delay - 1) + 1] would suffice. ■ 

The following theorem gives a sufficient field size for the 
required network error correcting (n, k) input convolutional 
code C s to be constructed with the required free distance 
condition (d/ ree (C s ) > 2t s + 1). 

Theorem 2: The code C s can be constructed and used to 
multicast k symbols to the set of sinks T along with the 
required error correction in the given acyclic, unit-delay, 
memory-free network with min-cut n (n > fc), if the field 
size q is such that 

n\(q - 1) 
and 

g >maJin 2m2[(ra+1)(T r-" 1) + 1] +2 



Proof: From the sufficient condition for the existence of 
a linear multicast network code for a single source network 
with a set of sinks T, we have 

q>\T\. 

Now we prove the other conditions. From the construction 
in [19], we know that a (n, fc, S) MDS convolutional code can 
be constructed over ¥ q if 



n\(q - 1) 



and 



q > 



dn 2 



Thus, with 5 = 2rk [(n + 1) (T de i ay ~ 1) + 1] as in Corol- 
lary Q] an input MDS convolutional code C s can be constructed 



over ¥ q if 



n\(q — 1) and q > 



2rn 2 [(n + I) (Tdeiay - I) + 1} 



+ 2. 



k (n — fc) 



n — k 

Such an MDS convolutional code 

the requirements in the construction 

(dfree(C s ) > 2m [(n + 1) (Tdeiay - 1) + 1] + 1 > 2t, + 1), 

and hence the theorem is proved. ■ 
2) Bound on Td free (C s ): Towards obtaining a bound on 
Td frcc (C s ), we first restate the following bound proved in [16]. 

Proposition 2: Let C be a (c, b, 5) convolutional code. Then 

T dfree (C) < (dfree(C)- 1)6+1. (13) 

Thus, for a network error correcting MDS convolutional 
code C s for the unit-delay network, we have the following 
bound on T dfrcc (C s ). 

Corollary 2: Let the CNECC C s be a (n, k,S = 
2rk [(n + 1) (T de iay - 1) + 1]) MDS convolutional code, 
where r and Tdeiay are as in Lemma [2] Then 

T dfTee (C*) < ^ nk [(« + 1) tdeiay - 1) + if 

+ 2rfc (n - fc) [(n + 1) (T d eiay - 1) + 1] + 1. 
Proof: For MDS convolutional codes, we have 

dfree(C) = (n-k)([S/k\ + 1) + 6 + 1 

With 5 = 2rk [(n + 1) (T d eia V - 1) + 1], we have 

dfree{C s ) = (fl - fc) {2r [(n + 1) (T de lay - 1) + 1] + 1} 

+ 2rfc [(n + 1) (T delay - 1) + 1] + 1 

dfree(C s ) = 2m [(n + 1) (Tdelay -l) + l]+7l — fc + l 

Substituting this value of d/ ree (Cs) and S in dot , we have 
proved that 

T dfrs ACs) < 4r 2 nk [(n + 1) (T delay - 1) + l] 2 

+ 2rfc (n - fc) \(n + 1) (T d eia V - 1) + 1] + 1. 



IV. Illustrative examples 
A. Code construction for a modified Butterfly network: 

Let us consider the modified butterfly network as shown 
in Fig. |2] with one of the edges at the bottleneck node (of 
the original unmodified butterfly network) having twice the 
delay as any other edge, thus forcing an inter-generation linear 
combination at the bottleneck node. The local kernels at the 
node defining the network code are the same as in that of 
the instantaneous butterfly case. We assume the network code 
to be over F2 and we design a convolutional code over F2 
that will correct all single edge errors in the network, i.e, all 
network error vectors of Hamming weight utmost 1. 

For this network, the matrix A is a 2 x 10 matrix having 
a 2 x 2 identity submatrix at the columns corresponding to 
edges ei and e%, and having zeros everywhere else. We assume 
B Tl and B T2 are 10 x 2 matrices such that they have a 2 x 2 
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and 




Fig. 2. Modified butterfly network with global kernels 

identity submatrix at rows ee, e% and eg, eirj respectively. With 
the given network code, we thus have the network transfer 
matrices at sink T\ and T2 as follows 

v ^3 



M Tl (z) 



where 

and 

where 

Ft 2 (z) -- 





z 4 z 2 z 3 z 2 



AF Tl (z) 



1 
z 1 



M T2 (z) 



AF T2 (z) 



z 








1 







z 
For single edge errors, we have the error pattern set to be 

$ = {{ ei }:i = 1,2, ...,9,10}. 

And thus the set W$ is the set of all vectors F2 that have 
Hamming weight utmost 1. The sets Wti and Wt 2 as in (fl4] i 
and $15[ at the top of the next page. Now 



and 



M Ti \z) 



M-Hz) 







z 
-A 



To obtain the processing matrices Pp 1 (z) and Pt 2 (z), let us 
choose the processing functions p T (z) = z A andp T (z) = z 3 . 



Then we have 

P Tl (z)=p Ti (z)M^(z) = 



z 




(16) 



Pt 2 (z)= Pt (z)M t 1(z) = 



1 



(17) 



j t 2 V'J^ ± T 2 

Therefore, W s can be computed to be as in ([Tol l at the top 
of the next page. Thus we have t s = 2, which means that we 
need a convolutional code with free distance at least 5. Let 
the chosen input convolutional code C s be generated by the 
generator matrix 

G I (z)=[l + z 2 1 + z + z 2 ]. 

This code has a free distance df ree (C s ) = 5 and Tdf ree (C s ) = 
6. Therefore this code can be used to correct single edge 
errors in the butterfly network as long as consecutive errors 
are separated by 6 network uses. With this code, the output 
convolutional code Ct x at sink T\ is generated by the matrix 

Go, Tl (z) = [z + z 3 



z 3 + z 4 + z 6 ] 



Now C Tl has d free (C Tl ) = 5 and T dfree (C Tl ) = 9 > 
2d /ree (C s ). As condition © is not satisfied, Case-B applies 
and hence the sink T\ has to use the processing matrix Pr 1 (z), 
and then decode on the trellis of the input convolutional code. 
Upon performing a similar analysis for sink T%, we have Table 
U as shown at the top of the next page. 

B. 4C2 combination network over ternary field 

We now give a code construction for double edge error 
correction in the 4C2 combination network with a network 
code over F3, shown in Fig. [3] with the given 2 dimensional 
network code, the network transfer matrices and the processing 
matrices (upon choosing the processing functions PTi(z) = 
Pt{z) = z V 1 < i < 6) corresponding to the 6 sinks are 
indicated in Table [II] 




Fig. 3. 4C2 unit-delay network 

The matrix Fr ± (z) corresponding to sink T\ is the 2 x 16 
matrix as follows 

M _ r z 1 1 ' 

Tl ^ Z) ~ z 1 
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W Tl = {(0, 0), (0, 1), (1, 0), (0, z), (0, z 2 ), (0, z 3 ), (0, z 4 ), (z, z 3 )} (14) 
W T2 = {(0,0),(0,l),(l,0),(z,0),(z 2 ,0),(z 3 ,0),(0,z 4 ),(z 4 ,z)} (15) 



W s = { (0, 0), (z 3 , z 2 ), (0, 1), (0, z), (0, z 2 ), (0, z 3 ), (0, z 4 ), (z, 0), (z 2 , 0), (z 3 , 0)} . (18) 



TABLE I 

Modified butterfly network with C s [df ree (C s ) = 5,T d/r(ji! (C s ) = 6] 



Sink 


Output convolutional code generator matrix [Go t; ( z )] 


dfree{CTi), ^d/reeC^O 


Decoding on 






z + z 3 z 3 + z 4 + z 6 




5,9 


Input trellis 


T 2 


[z 3 


+ z 4 + z B z + z 2 + z 3 ] 


6,12 


Input trellis 



For each sink, we have a similar Ft(z) matrix with a z-scaled 
identity submatrix and an identity submatrix and zeros at all 
other entries. 

For double edge error correction, the error pattern set $ is 



{{e;,ej} = 1, 2, 15, 16 and i ^ j} 



And therefore, we have the set W$ as the set of all 16 
length tuples from F3 with Hamming weight utmost 2. The 
set V i can be computed to be as shown in ( fT9l at the 
top of the next page. Now, the set 

W s ,t 2 = {w T2 {z)P T . 2 {z) I w T2 (z) G WrJ 

is computed to be as in (1201 1. also shown at the top of the next 
page. Similarly the sets W s .Ti (V 1 < i < 6) and 

W s = |J IV../ 

T;6T 

are computed. It is seen that for this network, 

= max wh (w s (z)) = 4 

■w s (z)ew s 

and 

= max (w T . (z)) =2, V 1 < i < 6. 

k> t . (z)GW Ti s 

Therefore we need a convolutional code with free distance 9 
to correct such errors. Let this input convolutional code C s 
over F3 be chosen as the code generated by 

Gi(z) =[l + z 2 + z 4 + z 5 2 + z + 2z 2 + 2z 4 + z 5 ] . 

This code is found to have df ree (C s ) = 9 with Td fr ^{C s ) = 
14. Thus it can correct all double edge network errors as long 
as consecutive network errors are separated by 14 network 
uses. The output convolutional codes T , c ; /ree (CT i ), their free 
distance and Td frcc (CTi) are computed and tabulated in Table 
HI1 at the top of the next page. For this example, all the sinks 
satisfy the conditions © and dS) for Case-A of the decoding 
and therefore decode on the trellises of the corresponding 
output convolutional codes. 



V. Comparison between CNECCs for 

INSTANTANEOUS AND UNIT-DELAY, MEMORY-FREE 
NETWORKS 

In the following discussion, we compare the CNECCs for 
a given instantaneous network constructed in [16] and the 
CNECCs of Subsection IIII-BI for the corresponding unit-delay, 
memory-free network. 

With the given acyclic graph Q(V,£), we will compare 
the maximum Hamming weight t s of any rt-tuple, over ¥ q [z] 
(w s (z) € W s , where W s is as in (|5]l) in the case of the unit- 
delay, memory-free network with the graph Q and over F g 
(w s £ W s where W s is as in ([TTl i) in the case of instantaneous 
network with the graph Q. 

Consider some w s (z) G W s such that 

w s {z) = wF T {z)P T {z) = wp T (z)F T (z)M T 1 (z) 

= [w St i(z),w St 2(z),... ,w s , n (z)] (21) 

where p T (z) and Pt{z) indicate the processing function and 
matrix chosen according to (fT2l) for some sink T G T, 
and w S: i(z) G F g [z]. We have Mt(z)\ z= i = Mt and also 
Ft(z)\ z =i = Ft, the network transfer matrix and the Ft 
matrix of the sink T in the instantaneous network. Now, by 
OTh . we have the n-length vector w S: i nst corresponding to the 
error vector w as 



wFtM? 1 



w s (z)\z 



where 



p T (z)\ z =i = 



PtO)U=i 

Det (M T (z))\ z=1 



g(z)\z=i 

by O. Now Det (M T (z)) \ z=1 = Det (M T ) ^ since M T 
is full rank. Also, g(z)\ z= \ ^ for the same reason. Therefore, 
p T (z)\ z= i ^ 0. Thus we have 



w H (w s j nst ) < w H (w s (z)) . 



(22) 



Therefore a CNECC for an instantaneous network may require 
a lesser free distance to correct networks errors matching one 
of the given set of patterns $, while the CNECC for the 
corresponding unit-delay, memory-free network may require 
a larger free distance to provide the same error correction 
according to the construction of Subsection IIII-BI 
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W S ,T 2 = < 





r feo), 


(0,z), 


(1,0), (0,1), 


(22,0), 


(0,22), 


(2,0), ' 




(0,2), 




(2,22), (22,2), 


(22,22), 


(2 + 1,0), 


(2 + 2,0), 


= < 


(22 + 1,0), 


(22 + 2,0), 


(2,1), (2,2), 


(22,1), 


(22,2), 


(M), 




(1,2^), 


(2,2), 


(2,22), (0,2+1), 


(0,2 + 2), 


(0,22 + 1), 


(0,22 + 2), 




(1,1), 


(1,2), 


(2,1), 




(2,2), 


(0,0) 




(z,2z), 


(0,2), 


(1,2), 


(0,1), 


(22,2), 


(0,22), 




(2,1), 


(0,2), 


M), 


{z,z), 


(22,22), 


(22,0), 




- 1,22 + 2), 


(2 + 2,22+1) 


, (22+1,2 + 2), 


(2z + 2,z + 


1), (2,22 + 1 


), (2,22 + 2) 


(22,2 + 1), 


(22,2 + 2), 




(1,22), 


(2,2), 


(2,22), 


(0,^ + 1), 


(0,2 + 2), 


(0,22 + 1), 


(0,22 + 2; 


(1,0), 


(1,1), 






(2,2), 




(2,0), 




(0,0) 



(19) 



(20) 



TABLE II 

4C2 COMBINATION NETWORK WITH C s [df ree (C s ) = 9, Tj.^ (C s ) = 14] 



Sink 


Network transfer matrix 


Processing matrix 


Output convolutional code 
gen. matrix [Go.Ti( z )\ 


dfreeiCTi ), 

T d 1rBe (C Ti ) 


Decoding on 


Ti 


M Tl (z) = 






P Tl (z) = 


(l N 
V 1 / 




[z + z :i + z 5 + z 6 2z + z 2 + 2z 3 + 2z 5 + z e >] 


5,9 


Output trellis 


T 2 


M T2 (z) = 


w 




Pt 2 (z) = 


? 1 J \ 
V 1 / 




[z + z 3 + z 5 + z 6 z 2 + 2z 6 ] 


6,11 


Output trellis 


T 3 


Mt 3 (z)=\ 


' z z s 
v 2z , 


) 


Pt 3 {z) = 


? 2 2 \ 
V 1 j 




[z + z 3 + z 5 + z 6 2z + 2z 2 + 2z 3 + 2z 5 ] 


6,11 


Output trellis 


T4. 


M Tl {z) = 


: : 




P Tl {z) = 


(I 2 \ 
\ 2 




[2z + z 2 + 2z 3 + 2z 5 + z 6 z 2 + 2z 6 } 


7,12 


Output trellis 


T 5 


M Ts (z) = I 


' Z S 

v z 2z / 


) 


Pt 5 (z) = 


12 2 \ 

K 2 J 




[2z + z 2 + 2z 3 + 2z 5 + z 6 2z + 2z 2 + 2z 3 + 2z 5 } 


9,14 


Output trellis 


T 6 


Mt 6 (z)= 1 


' 2 Z { 

K z 2z / 


) 


Pt 6 (z) = 


? 2 2 \ 
V 2 1 J 




[z 2 + 2z 6 2z + 2z 2 + 2z 3 + 2z s ] 


6,13 


Output trellis 



An example of this case is the code construction for double 
edge error correction for the 4C2 combination instantaneous 
network in [16] and for the 4C2 unit-delay network in this 
paper in Subsection IIV-BI It can be seen that while for the 
instantaneous network, the maximum Hamming weight of any 
w s £ W s is 2, the maximum Hamming weight of any w s (z) G 
W s in the unit-delay network is 4. Thus a code with free 
distance 5 suffices for the instantaneous network, while the 
code for the unit-delay network has to have a free distance 9 
to ensure the required error correction as per the construction 
in Subsection IIII-BI 

It is in general not easy to obtain the general conditions 
under which equality will hold in ( 1221 . as both the topology 
and the network code of the network influence the Hamming 
weight of any element in W s . For specific examples however, 
this can be checked. An example of this case is given in 
between the single edge-error correcting code construction for 
the butterfly network (over F2) for the instantaneous case in 
[16] (the additional intermediate node, headier) = V3 = 
tail(e§), does not matter for the instantaneous case), and 
for the unit-delay case in this paper in Subsection IIV-AI 
In both the cases, we have t s — 2, which means that an 
input convolutional code with free distance 5 is sufficient to 
correct all single edge network errors. However, as we see in 
Subsection IIV-AI processing matrices with memory elements 
need to be used at the sinks for the unit-delay case, while 
the processing matrix in the instantaneous case is just the 



M T matrix which does not require any memory elements 
to implement. 

VI. Simulation results 

A. A probabilistic error model 

We define a probabilistic error model for a unit delay 
network Q(V,£) by defining the probabilities of any set of 
i (i < |£|) edges of the network being in error at any given 
time instant as follows. Across time instants, we assume that 
the network errors are i.i.d. according to this distribution. 

Prob. (i network edges being in error) = p l (23) 
Prob. (no edges are in error) = q (24) 

where 1 < i < \E |, and p, q < 1 are real numbers indicating 
the probability of any single edge error in the network and 
the probability of no edges in error respectively, such that 

<z + Eli'iP J = i- 

B. Simulations on the modified butterfly network 

With the probability model as in (|23]l and d24]i with \£\ = 10 
for the modified butterfly network as in Fig. [2] we simulate the 
performance of 3 input convolutional codes implemented on 
this network with the sinks performing hard decision decoding 
on the trellis of the input convolutional code. In the following 
discussion we refer to sinks T% and T2 of Fig.|2]as Sink 1 and 
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Sink 2. The 3 input convolutional codes and the rationality 
behind choosing them are given as follows. 

• Code C\ is generated by the generator matrix 

G h (z) = [l + z 1], 

with df ree (Ci) = 3 and Td fm (C\) = 2. This code is 
chosen only to illustrate the error correcting capability of 
codes with low values of df ree (C) and Td free (C). 

• Code C2 is generated by the generator matrix 

G l2 {z)=[\ + z 2 1 + z + z 2 ], 

with df ree (C2) = 5 and Td fTCC (C2) = 6. This code cor- 
rects all double edge errors in the instantaneous version 
(with all edge delays being zero) of Fig. [2] as long as they 
are separated by 6 network uses. 

• Code C3 is generated by the generator matrix 

G h (z)=[l + z + z A l + z 2 + z 3 + z 4 ] , 

with d free (C 3 ) = 7 and T dfree (C 3 ) = 12. This code 
corrects all double edge errors in the unit-delay network 
given in Fig. [2] as long as they are separated by 12 
network uses. 

We note here that values of Td free (C) of the 3 codes 
are directly proportional to their free distances, i.e, the code 
with greater free distance has higher Td frss (C). Also we note 
that with each of these 3 codes as the input convolutional 
codes, the output convolutional codes violate at least one of 
the conditions of 'Case-A' of decoding, i.e, (Q,©, or (|9). 
Therefore, hard decision Viterbi decoding is performed on the 
trellis of the input convolutional code. 

Fig. [4] and Fig. [5] illustrate the BERs for different values 
for the parameter p (the probability of a single edge error) of 
( |23l l. Clearly the BER values fall with decreasing p. 

It may be observed that between any two of the 3 codes, 
say Ci and Cj — 1, 2, 3) there exist a particular value of 
p = pij where the BER performance corresponding to the two 
codes gets reversed, i.e, if code d has better BER performance 
than Cj for any p > then Cj performs better than Ci for 
any p < pij. Although such a cross-over value of p exists for 
each pair of codes, we see that all 3 codes have approximately 
the same crossover p value in Fig. |4](p ~ 0.16) and similarly 
in Fig.0(p« 0.15). 

With respect to such crossover points between the two codes 
Ci and Cj, we can divide the performance curve into two 
regions which we call as 'Td fTCC dominated region' (p values 
being greater than the crossover p value) and 'd/ ree dominated 
region' (p values being lesser than the crossover p value), 
indicating the parameter which controls the performance of 
the codes in each of those regions respectively. Again, because 
of the 3 crossover points being approximately equal to one 
another in each of Fig. [4] and Fig. [5] we divide the entire 
performance graph of all the 3 codes into two regions. The 
following discussion gives an intuition into why the param- 
eters Td frec {C) and df ree (C) control the performance in the 
corresponding regions. 

• df ree dominated region: In the df ree dominated region, 
codes with higher free distance perform better than those 



with less free distance. We recall from Proposition Q] 
that both the Hamming weight of error events and the 
separation between any two consecutive error events are 
important to correct them. Because of the fact p is low 
in the df ree dominated region, the Hamming weight of 
the modified error sequences of ( fTOb is less, and the 
error events that occur are also separated by sufficient 
number of network uses. Therefore the condition on the 
separation of error events according to Proposition Q] is 
automatically satisfied even for large Td }TCC (C) codes. 
Therefore codes which have more free distance (though 
having more T dsm (C)) correct more errors than codes 
with low free distance (though having less Td frcc (C)). It 
is noted that in this region the code C3 (which was de- 
signed for correcting double edge errors on the unit-delay 
network) performs better than C2 (which was designed 
for correcting double edge errors on the instantaneous 
version of the network). 

• Td free dominated region: In the Td free dominated region, 
codes with lower Td free (C) perform better than codes 
with higher T ( j /rcc (C), even though their free distances 
might actually indicate otherwise. This is because of 
the fact that the error events related to the modified 
error sequences of dTob occur more frequently with lesser 
separation of network uses (as p is higher). Therefore the 
codes with lower Td Sree (C) are able to correct more errors 
(even though the errors themselves must accumulate less 
Hamming weight to be corrected) than the codes with 
higher Td free (C) which demand more separation in net- 
work uses between error events for them to be corrected 
(despite having a greater flexibility in the Hamming 
weight accumulated by the correctable error events). 

Remark 2: The difference in the performance of code C\ 
between Sink 1 and Sink 2 is probably due to the unequal 
error protection to the two code symbols. When the code is 
'reversed' ,i.e. with Gi 1 (z) = [1 1 + z], it is observed 
that the performance at the sinks are also interchanged for 
unchanged error characteristics. 

VII. Concluding remarks 

In this work, we have extended the approach of [16] to 
introduce network error correction for acyclic, unit-delay, 
memory-free networks. A construction of CNECCs for acyclic, 
unit-delay, memory-free networks has been given, which cor- 
rects errors corresponding to a given set of patterns as long 
as consecutive errors are separated by a certain number of 
network uses. Bounds are derived on the field size required for 
the construction of a CNECC with the required error correction 
capability and also on the minimum separation in network 
uses between any two consecutive network errors. Simulations 
assuming a probabilistic error model on a modified butter- 
fly network indicate the implementability and performance 
tractability of such CNECCs. The following problems remain 
to be investigated. 

• Investigation of error correction bounds for network error 
correction in unit-delay, memory-free networks. 

• Joint design of the CNECC and network code. 
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• Investigation of distance bounds for CNECCs. 

> Design of appropriate processing matrices at the sinks 
to minimize the maximum Hamming weight of the error 
sequences. 

> Construction of CNECCs which are optimal in some 
sense. 

• Further analytical studies on the performance of CNECCs 
on unit-delay networks. 
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Appendix A 
Convolutional codes-Basic Results 

We review the basic concepts related to convolutional codes, 
used extensively throughout the rest of the paper. For q, power 
of a prime, let F g denote the finite field with q elements, ¥ q [z] 
denote the ring of univariate polynomials in z with coefficients 
from F q , ¥ q (z) denote the field of rational functions with 
variable z and coefficients from ¥ q and F g [[z]] denote the 
ring of formal power series with coefficients from ¥ q . Every 
element of F q [[z]] of the form x(z) — 2~2"iLo Xizt ' Xi ^ 



Thus, ¥ q [z] C F 9 [[z]]. We denote the set of n-tuples over 
F 9 [[z]] as F™[[z]]. Also, a rational function x(z) — |^|y with 
6(0) ^ is said to be realizable. A matrix populated entirely 
with realizable functions is called a realizable matrix. 

For a convolutional code, the information sequence u = 
[uq,Ui, ...,Ut] (ui G ¥ h q ) and the codeword sequence (output 



sequence)?; = [vo, Vi, v t ] (v, £ 
terms of the delay parameter z as 



can be represented in 



"0) = '■ 



and 



v(z) = ^2v. t 



i=0 



i=0 



Definition 4 ( [20]): A convolutional code, 
b/c (b < c) is defined as 



C of rate 



C = {v(z) e ¥ c q [[z}} | v(z) = u{z)G{z)} 

where G(z) is a b x c generator matrix with entries from ¥ q (z) 
and rank b over ¥ q (z), and v(z) being the codeword sequence 
arising from the information sequence, u(z) £ F£[[z]]. 

Two generator matrices are said to be equivalent if they 
encode the same convolutional code. A polynomial generator 
matrix [20] for a convolutional code C is a generator matrix for 
C with all its entries from F 9 [z] . It is known that every convo- 
lutional code has a polynomial generator matrix [20]. Also, a 
generator matrix for a convolutional code is catastrophic [20] 
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if there exists an information sequence with infinitely many 
non-zero components, that results in a codeword with only 
finitely many non-zero components. 

For a polynomial generator matrix G{z), let gij(z) be the 
element of G(z) in the i th row and the j th column, and 

Vi := max deg(gij (z)) 
j 

be the i th row degree of G(z). Let 

6 

S := v i 
i=i 

be the degree of G(z). 

Definition 5 ( [20] ): A polynomial generator matrix is 
called basic if it has a polynomial right inverse. It is called 
minimal if its degree 5 is minimum among all generator 
matrices of C. 

Forney in [21] showed that the ordered set {y\, v%, Ub} 
of row degrees (indices) is the same for all minimal basic 
generator matrices of C (which are all equivalent to one 
another). Therefore the ordered row degrees and the degree 
5 can be defined for a convolutional code C. A rate b/c 
convolutional code with degree 5 will henceforth be referred 
to as a (c, b, 5) code. Also, any minimal basic generator matrix 
for a convolutional code is non-catastrophic. 

Definition 6 ( [20] ): A convolutional encoder is a physical 
realization of a generator matrix by a linear sequential circuit. 
Two encoders are said to be equivalent encoders if they 
encode the same code. A minimal encoder is an encoder with 
the minimal number of delay elements among all equivalent 
encoders. 

The weight of a vector v(z) € F^[[z]] is the sum of the 
Hamming weights (over ¥ q ) of all its F^-coefficients. Then 
we have the following definitions. 

Definition 7 ( [20]): The free distance of a convolutional 
code C is given as 

df ree (C) = rain {wt(v(z))\v(z) G C, v(z) ^ 0} 



A method of constructing MDS convolutional codes based 
on the connection between quasi-cyclic codes and convolu- 
tional codes was given in [19]. The ordered Forney indices 
for such codes are of the form 

V\ = v 2 = ... = Vl < v l+x = ... = v k . 

where v\ = \5/k\ and Vk = \5 /k\ +1. 

It is known [19] that the field size q required for a (n, k, S) 
convolutional code C with df ree (C) meeting the generalized 
Singleton bound in the construction in [19] needs to be a prime 
power such that 

n 2 

n ^- 1 ^ M ^ S Kn-k) +2 - (25) 



Appendix B 

MDS CONVOLUTIONAL CODES 

We discuss some results on the existence and construction 
of Maximum Distance Separable (MDS) convolutional codes. 

The following bound on the free distance, and the existence 
of codes meeting the bound, called MDS convolutional codes, 
was proved in [18]. 

Theorem 3 ( [18]): For every base field F and every rate 
k/n convolutional code C of degree 6, the free distance is 
bounded as 

dfree(C) < (t» - k)([5/k\ + 1) + S + 1 

Theorem [3] is known as the generalized Singleton bound. 

Theorem 4 ( [18]): For any positive integers k < n, S and 
for any prime p there exists a field ¥ q of characteristic p, and 
a rate k/n convolutional code C of degree S over F g , whose 
free distance meets the generalized Singleton bound. 



